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Abstract. We present optical and infrared observations taken with the Very Large Telescope of the eclipsing binary pulsar 
system PSR B1718-19. The candidate companion of the pulsar, identified earlier in Hubble Space Telescope observations, has 
been detected in all three bands, R, I, and J. These detections allowed us to derive constraints on temperature, radius, and mass, 
pointing to a companion that has expanded to a radius between one of a main sequence star and one at the Roche-limit. We 
focus on the role of tidal circularisation in the system, which will have transformed the initially eccentric orbit expected from 
formation scenarios into the nearly circular orbit presently observed. Based on simple energy balance arguments, we are able 
to draw a picture of the companion's evolution resulting from the energy deposition in the star due to circularisation. In this 
picture, our measurement of the companion's parameters is consistent with the expected initial eccentricity. However, with the 
present understanding of tidal dissipation it remains difficult to account for the short time in which the system was circularised. 
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1. Introduction 

The binary radio pulsar PSR B1718-19 has a 1-s spin period 
and was discovered by Lyne et al. (1993) in the direction of 
globular cluster NGC 6342. The pulsar is young and has a high 
magnetic field: from the pulse period and its derivative one in- 
fers r char = P/(2P) = 10 Myr and B = 1.5 x 10 12 G. The binary 
orbit is circular, has a period of 6.2hr, and is accompanied by 
eclipses at low radio frequencies, of 400-600 MHz. The eclips- 
ing material is probably the companion's wind. 

The formation scenarios proposed (Lyne et al. 1993; Ergma 
1993; Wijers & Paczynski 1993; Zwitter 1993; Burderi & King 
1994; Ergma et al. 1996) fall into two main categories: (i) par- 
tial recycling of an old neutron star in a close encounter with 
other stars in the cluster core, and (ii) accretion-induced col- 
lapse (AIC) of a white dwarf. Van Kerkwijk et al. (2000) de- 
tected a candidate for the pulsar companion. The object's faint- 
ness was a sttong indication for scenario (i), though scenario 
(ii) could not be excluded definitively. 

Wijers & Paczynski (1993) noted that either formation sce- 
nario leads to an eccentric binary orbit, and tidal circularisa- 
tion must have worked to reduce the eccentricity to the cur- 
rently observed low value (e < 0.005). Thus, PSR B1718-19 
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might serve as an observational constraint on circularisation 
theory in a parameter region different from the region covered 
by constraints that are currently available: tidal cutoff periods 
of solar-type binaries (see, e.g., Mathieu et al. 2004 and ref- 
erences therein), binaries containing a giant star (Verbunt & 
Phinney 1995), or planets (Wu 2003). 

The circularisation efficiency depends strongly on the ob- 
ject's radius, but this is poorly constrained by the single-band 
detection of Van Kerkwijk et al. (2000). This motivated us to 
collect further observations to obtain a direct measurement of 
temperature and radius. In this paper, we describe the results. 
In §2, we describe the reduction and show our results in colour- 
colour and colour-magnitude diagrammes. We use these in §3 
to put further consttaints on mass, radius and temperature of 
the companion, and in §4 we investigate the implied role of 
tidal circularisation, focusing on the circularisation timescale 
and the consequences of tidal dissipation in the companion. 

2. Observations 

The observations on which we report here were done for us 
in Service Mode in June 2001 with the Very Large Telescope 
(VLT) at ESO (Paranal) through filters R, I and J (see Table 1). 
In this section we will describe their reduction and photometry. 

We will also use the data presented in Van Kerkwijk et al. 
(2000): images obtained through the F702W filter with the 
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Table 1. Log of observations 



Instrument 


Band 


Night 


tat (s) 


seeing (") 


#EXP 


STD-field" 


A sec z b 


VLT-FORS2 


I 


2001 Jun 20 


257 


0.48-0.61 


8 


PG1323 


0.04 






2001 Jun 22 


240, 257 


0.44-0.67 


8 


SA92-249 


0.00 


VLT-FORS2 


R 


2001 Jun 20 


317 


0.51-0.68 /0.75 c 


11/12" 


PG1323 


0.02 






2001 Jun 22 


300, 317 


0.48-0.66/ 0.7 l-0.83 c 


11 / 13 C 


SA92-249 


0.02 


VLT-ISAAC 


J 


2001 Jun 7 


30.30 


. . . / 0.59-0.87 r 


0/W 


S860-D 


0.00 






2001 Jun 14 


30.30 


0.36-0.56 


48 


S860-D 


0.07 


H S T-WFPC2 


F702W 


1997 Mar 8 


700 


0.046 


12 







" The standard magnitudes are taken from Stetson (2000) for I and R and from Persson et al. (1998) for J. 

* This column lists the difference in airmass between the single exposure science frame that was used for calibration and the standard frame 
(see text). 

c The frames were selected on seeing value. Listed in the seeing column is, before the slash, the seeing range of frames included in the 
mosaic image, and, after the slash, the seeing range of frames that were excluded. In the # EXP column the number of exposures included in 
the mosaic image is displayed, with after the slash the total number of exposures taken. 



Wide Field and Planetary Camera 2 (WFPC2) of the Hubble 
Space Telescope (HST) in March 1997. Only the data from the 
Planetary Camera (800x800 pixels, at 07046 pixel -1 ) will be 
used for the present paper. As will be discussed in §2.3, we 
were unable to derive reliable photometry from these data. The 
F702W images still serve, however, for measurements of the 
proper motions of stars in the field. 

The astrometry done by Van Kerkwijk et al. (2000) enables 
us directly to identify the candidate for the pulsar's companion 
in our images, shown in Fig. 1 . 

2.1 . Infrared observations in J 

The J-band observations were done with the short-wavelength 
channel of the Infrared Spectrometer And Array Camera 
(ISAAC; Moorwood et al. 1998), mounted on the 8.2m Antu 
(Unit Telescope 1). The detector was a Rockwell Hawaii 
1024x1024 pixels array, with a plate scale of 07148 pixel -1 . 

The ISAAC-observations need to be corrected for two in- 
strumental effects: the odd-even column effect, which can be 
seen as an offset between the odd and even columns of the ar- 
ray, and the electrical ghost, which consists of an additional 
signal proportional to the sum of the intensity along a given row 
and the row 512 rows away (see the ISAAC Data Reduction 
Guide 1 , Amico et al. 2001). We employ ESO's software pack- 
age eclipse (Devillard 1997) for these corrections. Next, we use 
midas to subtract the dark current from all frames and divide 
them by twilight flat fields. 

From the images, it was clear that the seeing in the first 
observing night was much worse than in the second night (see 
Table 1). We decided to exclude all 16 frames taken on the 
first night. We register the 48 frames of the second night with 
respect to each other, applying offsets rounded to the nearest 
integer pixel, and compute an average image. Bad pixels are 
excluded from this average. We regard a pixel as bad, if (i) it is 
located at rows 513 or 514 (2048 pixels), where bias-level vari- 
ations with time prevent proper bias subtraction; (ii) the differ- 

1 Available at http://www.eso.org/instruments/isaac/drg/drg.ps 



ence between a long and a short dark-current exposure at that 
location lies more than 7.5cr away from the mean (57 pixels); 
(iii) it has a value lower than 2.75<x below background in more 
than 50% of the frames (172 pixels), or (iv) if it is a hot pixel 
(59 pixels). Finally, in order to prevent problems with differ- 
ences in the exposure time and distortion near the edges, we 
extract the central 600 x 600 pixels (89 x 89") part of our mo- 
saic image. In this subimage, the point-spread function (PSF) 
shows negligible variation. 

We perform PSF-fitting photometry using the DAOPHOT 
II package (Stetson 1987). The determination of the PSF is 
hampered by the crowdedness of the field: one object per 4 
arcsec 2 , at a full width at half maximum (FWHM) of 075. We 
deal with this following the suggestions of Stetson (1987), by 
first estimating the PSF based on several relatively isolated 
stars, and then using this PSF to subtract neighbouring stars. 
On the new image a better PSF can be determined. This proce- 
dure has been repeated four times. With the resulting final PSF, 
the magnitudes of all stars are determined. 

Observations of a standard star (Table 1) serve to calibrate 
the photometry. As these observations are taken with very short 
exposure times, the PSF varies greatly over the field. Hence, 
PSF-fitting is not reliable and we resort to aperture photome- 
try. The standard star's aperture magnitudes can then be related 
to magnitudes of stars in a single exposure science frame, mea- 
sured with a corresponding aperture. For this purpose, we select 
a science frame taken at airmass as close as possible to that at 
which the standard star was taken (see Table 1), so that there 
is no need to correct for atmospheric extinction. To account for 
the quite different seeing of the standard and the science ex- 
posures, we use apertures with an equal radius in units of the 
image FWHM, not pixels. We used the trend of magnitude as 
a function of aperture to verify that using an aperture of the 
same radius in FWHM indeed results in a superior match of 
the fraction of the flux inside a given aperture in two different 
frames. 

Next, we compute an aperture correction to relate the aper- 
ture magnitude in the single science frame to the magnitudes 
determined from PSF fitting in the mosaic image. Out of ne- 
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Fig. l.Images of the field of PSR 
B 1718—19. Shown are - clockwise - 
the images in R and I (taken with VLT- 
FORS2), in F702W (ffiST-PC), and in J 
(VLT-ISAAC). The F702W-observations 
have been presented in Van Kerkwijk et al. 
(2000), who also performed the astrometry 
to derive the 95% confidence error circle 
of radius O".^, centered on the position of 
the pulsar. This error circle is displayed 
in the images, and contains a detection 
of the candidate counterpart in all bands. 
For stars indicated by numbers in the 
F702W-panel, and for the PSR-companion 
we list positions and magnitudes in Table 2. 



Table 2. Positions and apparent magnitudes of the candidate compan- 
ion of PSR B1718-19 and surrounding stars. The uncertainty in the 
photometric zero point is 0.02 mag for J and < 0.03 mag in R & I. 
The uncertainty in the astrometric tie is 0'.'06 in each coordinate (Van 
Kerkwijk et al. 2000). Photometry for all stars in the field, including 
uncertainties, is available online. 



ID a 2 ooo <?2000 R I J 

17 h 21 m -19°36' 



PSR 


01 s :548 


36'.' 77 


25.87 


24.25 


22.59 




±0:004 


±0706 


±0.20 


±0.08 


±0.08 


1 


01 s :667 


37/55 


23.31 


22.09 


20.56 


2 


01 s .421 


37/26 


24.79 


22.98 


21.08 


3 


01 s :625 


32'.' 91 


21.00 


20.35 


19.37 


4 


01 s :608 


41'.' 54 


25.74 


23.96 


22.03 


5 


01 s :828 


32'.' 87 


23.55 


22.53 


21.14 


6 


or.m 


35/67 


20.50 


19.82 


18.75 


7 


01U56 


34'.' 34 


24.82 


23.49 


21.67 


8 


oi s : 114 


32'.' 06 


20.02 


19.40 


18.49 


9 


01 s :767 


44'.' 39 


19.59 


18.95 


17.99 



2.2. Optical observations in R and I 

The R and I-band observations were carried out with the FOcal 
Reducer/low dispersion Spectrograph (FORS2; Seifert et al. 
2000) mounted on Yepun (VLT-UT4). FORS was used in 
Standard Resolution imaging mode, in which the plate scale 
on the 2048 x 2048 pixels SITe detector is 07201 pixeF 1 . 

The reduction was done in a similar fashion as described 
above for the ISAAC-observations. First, all frames are bias- 
subtracted and divided by twilight flat fields. We select frames 
with a seeing lower than 077 (see Table 1) and compute an av- 
erage mosaic of the registered frames, excluding 1329 pixels in 
each image located on a bad column, as well as pixels that were 
hit by cosmic rays. From the mosaics, we extract the central 
400 x 400 pixels (80 x 80") parts, which are covered by all in- 
dividual frames and which have negligible distortion. On these 
subimages the photometry is done using procedures described 
above for the ISAAC data. For the calibration, we determine a 
weighted average from both observing nights for the zero point 
and the color term. This leads to a total error in the calibration 
of ^0.03 mag in both bands. 



cessity, since we do not have another infrared colour, we ne- 
glect any colour terms. From the ISAAC Data Reduction Guide 
(Amico et al. 2001, §5.5.1), this appears reasonable: the J - K 
colour term is expected to differ from by less than 0.01. 
Overall, we estimate that the total error in the calibration is 
~ 0.02 mag. 



2.3. Problems with the F702W-observations 

We tried to obtain photometry from the HST F702W observa- 
tions, both using HSTphot (Dolphin 2000) and following the 
instructions in the HSTWFPC2 data handbook (Baggett 2002). 
We failed to obtain numbers that we felt were sufficiently reli- 
able, for two reasons. First, indepent of the reduction method, 
we found that in colour-colour diagrammes (e.g., m F7 02w - I 
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vs. I - J) the root-mean-square scatter for bright, blue stars 
(I-J < 1.05) was as large as 0.026 mag, while in a comparable 
plot for R-I (see Fig. 2), it is only 0.008 mag. This implies the 
WFPC2 data had a base noise of ~0.02 to 0.03 mag, which we 
do not understand. From discussions with Drs Andrew Dolphin 
and Peter Stetson (2004, personal comm.), we learned that 
such scatter is normally present, but that it should average out 
when one uses multiple frames. In our case, however, it did 
not, possibly because our dither pattern used integer pixel off- 
sets. We also confirm it is present in other comparisons be- 
tween ground-based and WFPC2 data (such as those presented 
in van Kerkwijk & Kulkarni (2001); for that field, we find that 
ACS/HRC data show much less scatter and are completely con- 
sistent with the ground-based data). 

The second and more severe problem is that the two re- 
duction schemes lead to a zero point that differs by 0.1 mag. 
We are confident that this is not because we did an inappro- 
priate aperture correction (such as the 0. 1 mag offset between 
the standard aperture of 075 and "nominal infinity"); we do not 
find such differences for images in other filters that we anal- 
ysed. Since we were unable to identify the flaw in the reduction, 
we conclude that we cannot rely on the magnitude derived for 
the pulsar companion, and we excluded the F702W photometry 
from our analysis. We note, however, that within a 0.1 mag un- 
certainty, the magnitude found is consistent with that expected 
from our R and I-band detections. 

2.4. Photometry and Proper Motions 

We used our calibrated photometry (Table 2 and on-line ma- 
terial), to construct colour-colour and colour-magnitude dia- 
grammes; see Fig. 2. In the colour-magnitude diagramme, we 
find a rather broad main sequence, which suggests that many 
of the stars do not belong to the cluster. 

As the HST observations were taken 4.25 years prior to 
the VLT-observations, we can try to separate cluster stars from 
fore- and background objects using proper motions. For this 
purpose, we measured positions in the HST PC and our best- 
seeing VLT observation, that in the J band. The result is shown 
in Fig. 3. One sees that especially among brighter stars, a clus- 
tering of stars associated with NGC 6342 is apparent. This is 
confirmed in Fig. 2, where we use filled symbols to indicate 
stars that are likely cluster members (those within the dashed 
circle in Fig. 3): the brighter ones, with J < 21, form a clear, 
narrow sequence that is well matched by the expected main se- 
quence for NGC 6342. 

Nearly all bright stars that have proper motions inconsis- 
tent with cluster membership (outside the dashed circle; open 
symbols in Fig. 2) are brighter at a given colour than the cluster- 
associated ones. This is expected, since at the galactic latitude 
of 9°.73 of our observations, most stars belonging to the disk 
should be in front of the cluster. We can estimate a typical dis- 
tance by considering the number of stars per square arcsecond 
at a given distance. For an exponential distribution, this should 
scale as n(d) oc d 2 exp(-dsin\b\/zo), which reaches a maximum 
at c/ max = 2zo/sin\b\. Given a latitude b = 9°J3 for NGC 6342 
(Harris 1996, February 2003 revision) and scale height zo = 




1 1.5 
I - J 




12 3 4 

R - J 



Fig. 2. a Colour-colour and b colour-magnitude diagrammes of the 
companion of PSR B 17 18— 19 (error bars) and stars in the field. The 
filled circles denote stars with a proper motion inside the dashed circle 
of Fig. 3, which are likely associated with the cluster NGC 6342. The 
open circles indicate stars with a proper motion outside this circle. For 
stars denoted by a dot, no proper motion is available, since they were 
outside the field of view of the PC. The detection limit in J is around 
23 mag. The solid line represents model main sequence stars, labeled 
with their masses (in solar units), and - in italic typeface - temperature 
(small tickmarks denote intervals of 100 K). The models are corrected 
for a reddening £ v _/=0.55 and distance modulus (m-M) = 14.69 (see 
§3). The dotted arrow denotes a reddening of E V -i =0.55 ± 0.07 (with 
the size of the solid part equal to the uncertainty). 



325 pc (Reid & Majewski 1993), we obtain d mix - 3.8 kpc. 
Since the distance of NGC 6342 is 8.6 kpc, we would thus ex- 
pect stars in the galactic disks to be typically ~ 1 .7 mag brighter 
than stars in NGC 6342. This is roughly consistent with the off- 
set seen in Fig. 2b. 

From Fig. 3, one sees that for stars with J > 21, the distinc- 
tion between cluster-associated and foreground stars is much 
less clear. This is due to the increasing uncertainty in the J-band 
position (which dominates the proper motion uncertainty). For 
stars with J < 20, the position's uncertainty is ~ 2 mas, while 
around J « 21, it has increased to ~ 4 mas, reaching up to 
~ 10 mas for stars at the detection limit (/ w 23). The latter er- 
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Fig. 3. Proper motion of the counterpart of PSR B 17 1 8- 1 9 (error bar) 
and stars in the field. The three different sizes of the dots correspond 
to three brightness classes,: J < 20.3, 20.3 < J < 21.7, and J > 
21.7. The dashed circle, with a radius of 2.3 mas yr -1 , indicates the 
group of stars which we associate with the cluster. The counterpart 
of PSR B1718-19 has a proper motion consistent with our estimate 
of the cluster's proper motion (the uncertainty in declination is much 
larger than that in right ascension due to the presence of a bright star 
to the south of the companion). 

ror corresponds to an error in the proper motion of 2.3 mas yr~ 1 , 
equal to the radius of the dashed circle in Fig. 3. 

3. Analysis 

In this section, we estimate the companion's radius and tem- 
perature from its observed magnitudes. Since it is associated 
with NGC 6342, we can use the distance modulus (m - M)o= 
14.69 ± 0.20 found by Heitsch & Richtler (1999). Heitsch & 
Richtler also determined the reddening, but found a rather large 
differential variation. Below, therefore, we first determine the 
reddening in our field, by fitting models to cluster stars in the 
colour-colour and colour-magnitude diagrammes. 

Both to determine the reddening and to infer the compan- 
ion's properties, we use the evolutionary tracks for low-mass 
stars from Baraffe et al. (1998) (hereafter, BCAH models). We 
use the tracks for metallicity [M/H] = -0.5, which corre- 
sponds to [Fe/H] -0.77 (Baraffe et al. 1997), close to the 
measured metallicity [Fe/H] = -0.65 for NGC 6342 (Harris 
1996, February 2003 revision). 

3.1. Reddening 

Heitsch & Richtler (1999) found strong evidence for differen- 
tial reddening for NGC 6342, with Ey-i varying between 0.64 
and 0.74 in the core region of the cluster. At the location of 
PSRB1718-19, they inferred E v _, = 0.70... 0.76. 



In Fig. 2a, we see that the colours of the likely cluster mem- 
bers are at best marginally consistent with such high reddening; 
a much better fit is found for a value E V -i — 0.55. To try to 
understand this discrepancy, we tried to determine whether the 
reddening varied over our smaller, 80 x 80" field. We find some 
evidence for variation, with Ey-i - 0.40 at the south end, and 
E V -i - 0.75 at the north end. Thus, it may be that the value 
found by Heitsch & Richtler (1999) simply reflects a different 
effective centre of the region in which the reddening was deter- 
mined. 

Looking in detail at Fig. 2, one sees that for masses above 
~ 0.6 M Q , the models are slightly too red in R - I. This could 
be corrected by decreasing the reddening further, but then the 
break in the curve at ~0.6 M Q would happen at an / - J colour 
that is too blue. The fact that we cannot obtain a better overall 
fit is likely related to inaccuracies in the models (see Bertone 
et al. 2004). To take this uncertainty into account, we adopt a 
relatively large uncertainty in Ey-i, of 0.07 mag. 

We conclude that the reddening to PSR B 17 18- 19 is 
E v -i = 0.55 + 0.07 (or E B . V = 0.40 + 0.05). Converting red- 
dening to extinctions in the various bands (Schlegel et al. 1998, 
Table 6), we find A R = 1.07, A/ = 0.78, and Aj = 0.36. 

3.2. Radius & Temperature 

Below we will determine the radius and temperature through a 
direct fit to the models. Before doing so, we will make a cruder 
but more insightful estimate. We looked at the BCAH-models 
and found that, to satisfying approximation, colors depend on 
temperature and not on surface gravity, as long as A log g < 1 . 

From Fig. 2a, we see that the companion's colours are con- 
sistent with those of main sequence stars with a temperature 
of ~ 3400 + 150 K. The influence of the error in the red- 
dening can be estimated by moving the star over the solid 
part of the reddening vector; this induces an error of ~ 70 K. 
Adding this in quadrature, the final estimate of the temperature 
is 3400 + 170 K. 

The companion's radius can now be estimated from Fig. 2b. 
The masses of main-sequence stars with temperatures consis- 
tent with that of the companion lie in the range 0.13 < M < 
0.27M©. Figure 2b shows that these stars are fainter than the 
companion, which implies that the companion has a larger ra- 
dius, by a factor lQ-°- 2AJ (where AJ is the brightness differ- 
ence). From Fig. 2b, we estimate AJ 0.3 mag and AJ 2.0 
mag for MS-stars of mass 0.27 and 0.13 M e , respectively. 
Given the main-sequence radii from the models (stars with 
mass (0.13, 0.25, 0.30) M have radii (0.15, 0.26, 0.30) R Q ) one 
infers a radius of ~ 0.3 1 R Q at the high temperature end, and 
0.37 R G at the low end. 

For a given temperature, the uncertainty in the radius thus 
derived has a contribution from the uncertainties in the observa- 
tions and the reddening (about 0.09 mag for the J band), but it is 
dominated by the 0.2 mag uncertainty in the distance modulus. 
The uncertainty in the radius for given temperature is therefore 
-10%. 

Now that we know what temperature, radius, and associ- 
ated uncertainties we expect from our observations, we pro- 
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ceed with a fit. We want to take into account the uncertainty in 
distance and reddening and hence define x 2 as 

_ / /!)* — A"" — An) 4- *i no(K H \ — \" 



-2- 

Y V 



4var . 



Am var + 5 log (flvar/Amod) 



M mod 



' pvar 



c y -i 



o~ EM 



o-x 



Am var - Amu, 



1 



(1) 



with A™ r 



O" Am.lit 

(Ax/^y-/)^^, and Am = (m - M) . The la- 
bel "var" indicates that that part of the equation can be var- 
ied to minimize x 2 , and the symbol X denotes the filter: R, I, 
or J. Thus, we have five measurements (three magnitudes rtiR, 
mi, and m j, the distance modulus Am\ lU and reddening £'"_,) 
and four parameters (radius and temperature [which enter the 
model magnitudes], true distance modulus Am var and reddening 
E™j). Using Eq. 1, we create a contourplot of x 2 as a function 
of radius and temperature in Fig. 4. Here, we have marginal- 
ized over reddening and distance modules (i.e., for every radius 
and temperature, reddening and distance modules are adjusted 
to minimize x 2 )- We see that the result is consistent with the 
crude estimate made above. 

From Fig. 4, we see that the companion likely has a ra- 
dius larger than that of a main-sequence star of the same mass. 
The companion mass is constrained by the pulsar mass func- 
tion (M c > 0.11 M Q for a 1.35 M pulsar, Lyne et al. 1993). 
Furthermore, for given mass, the radius of the companion has to 
be smaller than the Roche lobe, and the temperature cannot be 
cooler than the Hayashi limit, where a star is completely con- 
vective. The latter limit implies that for given mass, the com- 
panion should be on or hotward of the corresponding pre-main- 
sequence track. 

The range of allowed parameters of the candidate compan- 
ion is listed in Table 3. At the l-cr confidence level, we observe 
a star with 0.11 < M < 0.27 M that is bloated compared to 
a main-sequence star of the same mass. Its radius lies below 
the Roche-limit for the mass range M > 0.14 M Q . However, 
at the 2-cr level we cannot exclude either a somewhat more 
massive main-sequence star or a lower-mass, Roche-lobe fill- 
ing star (see Fig. 4). As we will see below, however, our most 
likely intermediate position may be most consistent with ex- 
pectations from the preceding history. 

4. Discussion 

We now turn to the implications of our measurements. In do- 
ing this, we will first assume that the system formed through 
the scenario of a three-body interaction in the core of the clus- 
ter. The other, in our opinion less likely scenario of accretion- 
induced collapse of a white dwarf will be discussed seperately 
in §4.3. 

If the system formed in a three-body interaction, a bloated 
companion is expected. This is because the binary that is pro- 
duced is expected to be highly eccentric, with e 0.7 (Phinney 
1992, §5), while the system's current eccentricity is very low 
(e < 0.005). Thus, the orbit must have circularised within the 
system age, which is 10 Myr at most (Van Kerkwijk et al. 
2000). The natural mechanism - tidal circularisation - trans- 
fers orbital energy to the companion, thereby bloating the star. 



0.5 



0.1 

4000 





i 


. — i 

-O 03 

°d ° 
. ■ . o 








•' / / 


















oTB^^s^k^^^X^^ 




















Wo / 






















32 


i 


o.sNl 

0.175N^ 
0.15 

i 


0.13 


< -. 40 
51 

-3^; 64 
•. 81 
102. 

0.11 



3800 



3600 3400 
Temperature (K) 



3200 



3000 



Fig. 4. Constraints on the radius and temperature of the companion of 
PSR B1718-19. The ellipses show the 1, 2, and 3-cr contour levels 
(Ax 2 = 2.30, 6.17, and 11.8). The bold continuous line connect- 
ing the filled circles labeled with stellar masses (in solar units) in- 
dicates the model main sequence stars (the BCAH-models with age 
10 Gyrs). The dotted lines connect BCAH-models in a pre-main- 
sequence (PMS) contraction phase; the bold dotted line marks the 
lower limit set to the companion mass by the pulsar mass function. 
These PMS-tracks are close to the Hayashi-limit, and stars in hydro- 
static equilibrium can be only on or to the left of them. The bold 
dashed line connects the Roche-radii for stars on the PMS-tracks. The 
hatched region shows the full range of allowed parameters in the 2cr 
confidence region (see Table 3). Finally, the thin continuous lines 
are isochrones labeled with their age in Myrs. These isochrones show 
that, within the system age of 10 Myrs, the companion cannot have 
shrunk back to the main sequence if any significant bloating occurred 
that involved the whole star. 



Below, we first estimate the expected bloating, and then 
look into whether tidal theory can account for the observed 
limit on the circularisation timescale. 



4.1. Energy balance during circularisation 

The response of a star (or planet) to the dissipation of tidal en- 
ergy has been investigated by a number of authors (e.g. Ray 
et al. 1987; Podsiadlowski 1996; Gu et al. 2004; Ivanov & 
Papaloizou 2004), but none of these studies is directly applica- 
ble to our situation. Specifically for PSR B 17 18- 19, Verbunt 
(1994) showed that the tidal energy is of the order of the poten- 
tial energy of the star, assuming a ini » a c or e ini ~ 1 (here and 
below, the labels "ini" and "c" refer to the situations where the 
pulsar binary has just been formed and where the orbit has been 
circularised, respectively). With our constraints on the compan- 
ion's radius and temperature, we can make a more precise pic- 
ture of the exchange of energy between orbit and stellar inte- 
rior. 
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Table 3. Companion parameters, at 2<x confidence level and at maximum probability, and corresponding circularisation timescales. 
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Position in 


Temperature 


Radius 


Mass a 


Tcirc [Z] 


r circ [G&Of 


R, T-diagram 


(K) 


(Re) 


(Mo) 


(Myr) 


(Myr) 


Main sequence 


3560 - 3700 


0.26 - 0.33 


0.25 - 0.35 


17-5.3 


(10- 3.4) x 10 3 


Maximum probability 


3400 


0.34 


0.15 


0.41 


230 


Roche-lobe 


3310- 3440 


0.38-0.43 


0.11-0.17 


0.088 - 0.089 


46-50 



" The mass of the companion is estimated in Fig. 4 through interpolation of the BCAH-models. 

* The circularisation timescales are computed in two prescriptions. One is described in Zahn (1989, indicated with [Z]), the other in Goodman 
& Oh (1997, indicated with [G&O]). 



For our estimate, we assume that at the onset of circulari- 
sation, the companion is an ordinary main sequence star. Given 
the range of masses allowed for the companion (see Table 3 
and Fig. 4), it will be fully convective. A measure for the 
timescale on which energy is distributed across the star is the 
convective friction time tf = (mR 2 /L) 1 ^ 3 , which varies from 
0.55 yr at the high mass end of allowed companion parame- 
ters to 0.43 yr at the low mass end. This is much shorter than 
the thermal timescale in most of the star and than the circu- 
larisation timescale. We thus assume that the energy dissipated 
during circularisation is distributed through the entire star. 

In reaction to the increase in internal energy, the star will 
expand until hydrostatic and approximate thermal equilibrium 
is regained. For a completely convective star, the total energy 
E+ is given by, 



1 



^pot 



3Gm 2 
1R 



(2) 



where we used the Virial theorem and assumed the star's po- 
tential energy £ pot was roughly that of a n = 1.5 polytrope 
(Chandrasekhar 1939, Eq. 90), as appropriate for a completely 
convective star. The symbols m and R denote the companion's 
mass and radius, respectively. 

Next, we use that orbital angular momentum should be 
roughly conserved in the circularisation process (the compan- 
ion's rotational angular momentum is much smaller than the 
orbital angular momentum). This implies a; n j(l - e? ; ) = a c (l - 
e 2 ) « a c , and hence, 



-^orbit.ii 



GmM 



£orbit,c(l 



eld, 



(3) 



where M is the mass of the neutron star, for which we take 
1.35M (Thorsett & Chakrabarty 1999), and a the semi-major 
axis, which can be computed for given companion mass us- 
ing the observed period of 22314.83 s (Lyne et al. 1993). 
Combining Eqs. 2 and 3, we can write the exchange of orbital 
energy and total energy of the companion, A£ or t,i t = A£*, as: 



GmM 
— 1 

2a c 



3Gm 2 



1 

RfAS 



(4) 



where we used the assumption that the initial radius R ini is that 
of a main-sequence star, R M $- Here, we neglected the contribu- 
tion from the companion's rotational energy, since this is much 
smaller than the change in orbital energy. 

For a given companion mass, we can use Eq. 4 to calcu- 
late the initial eccentricity that induces bloating of the main 



sequence star up to a given radius. In this way, we convert 
the lcr and 2cr contours in the radius-temperature diagramme 
(Fig. 4) into corresponding contours in an eccentricity-mass di- 
agramme (Fig. 5), where we determine the mass by assuming 
the star is presently on a cooling track. 

From Eq. 4 we can also calculate the critical initial ec- 
centricity, which is the eccentricity for which the main se- 
quence star is bloated exactly to the Roche limit. For a rough 
estimate, we approximate the semi-major axis, which is only 
weakly dependent on m, by a c « 2R Q , and the Roche-limit by 
R L * 0A6a(m/(m + M)) 1/3 (Paczyriski 1971). Applying the 
mass-radius relation R M s/R Q « m M s/M , we obtain 



1.270-0.553 



Q.2M, 



2/3 



(5) 



In Fig. 5, one sees that the approximation agrees rather well 
with the values determined directly from the model; the differ- 
ences are largely due to the inaccuracy in the assumed mass- 
radius relation. 

In principle, the companion could fill its Roche lobe not 
only at the end of circularisation, when the companion's radius 
is maximal, but also at earlier stages when the radius is smaller 
but the orbital separation at periastron is smaller as well. To find 
the radius at any value of e, we solve the differential equation, 



dE+(R) dE mbil (e) 



de 



de 



3Gm 2 d(l/R) _ GmM 
7 de~ 



«,(6) 



where again assumed conservation of orbital angular momen- 
tum. We find 



R(e) = 



Rms 



1 + (7MR MS )/(6ma c )(e 2 - e 2 ) 



(7) 



Requiring that this radius is smaller than the Roche limit (using 
the approximation of Eggleton 1983) throughout the circular- 
isation process, we can compute an upper limit to the initial 
eccentricity above which Roche lobe overflow occurs. This up- 
per limit is also shown in Fig. 5. 

Above, we neglected any energy loss during or after cir- 
cularisation. From the isochrones in Fig. 4 we conclude that, 
within the system age of 10 Myrs, the star cannot have 
shrunk back significantly due to radiative losses. A contraction 
timescale (E+/L) substantially longer than the tidal bloating 
timescale (E^/L^fe) is indeed expected, since the star's lumi- 
nosity of ~ 0.014 L is substantially smaller than the tidal lu- 
minosity L t ide ^ (e 2 ni GmM/2a c )/T = 0.08 L Q . This leaves only 
one other potential flaw in our approach: whether the energy 
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Mass (M Q ) 

Fig. 5. Initial eccentricity ej n i for which the change in orbital energy 
A£ orbil during circularisation equals the excess energy fi^obs ~ £*,ms 
of the present companion relative to a main-sequence star (see Eq. 4). 
Shown are the initial eccentricities that can account for the bloating 
observed at the best-fit point as well as at the 1 and 2-<x contours from 
Fig. 4. The hatched region indicates the lower limit to the companion 
mass of 0.11 M Q . The dotted line shows the estimate from Eq. 5 for 
the critical initial eccentricity at which the star would be bloated to the 
Roche lobe. The crosses mark the critical eccentricities found from 
a more careful calculation; the main difference is due to the use of 
the more accurate main-sequence radii from the BCAH-models. The 
bold dashed line indicates the critical eccentricity above which Roche 
lobe overflow occurs at any time during the circularisation process 
(see text). 



is dissipated all through the star. The dissipation could be con- 
fined to a thin layer (see the discussion in §5.2 of Van Kerkwijk 
et al. 2000), and turn the energy-transport beneath the layer ra- 
diative. This would require a much more sophisticated treat- 
ment, if only because of the implications for the dissipation 
mechanism itself. 

In summary, if our assumption that the dissipation energy is 
distributed throughout the whole star holds, Figure 5 shows that 
our measurement of the companion's temperature and radius is 
in good agreement with the initial eccentricity of ~0.7 expected 
from the formation model. 

4.2. Circularisation timescale 

Regardless of the dissipation mechanism, tidal circularisation 
is expected to produce a decay of the eccentricity that is ex- 
ponential. To reduce the eccentricity from the expected initial 
eccentricity of ~ 0.7 to the tiny currently observed value of ^ 
0.005 therefore requires at least four circularisation timescales. 
Since the system is likely to be in its present state for at most 10 



Myr (Van Kerkwijk et al. 2000), the circularisation timescale 
must have been ^ 2.5 Myr. 

A star of a mass in the observed range is expected to be 
fully convective. The dominant dissipation mechanism in a 
convective star is thought to be the friction between the con- 
vective motions and the tidal flow, produced by turbulent vis- 
cosity (see, e.g., Zahn 1989). This mechanism becomes less ef- 
ficient if - like in our case - the convective turnover timescale 
is longer than the tidal period, but it is not clear to what ex- 
tent. We use two prescriptions for the effect of the timescale 
mismatch: Zahn (1989), in which the efficiency is assumed to 
decrease linearly with the ratio convective-to-orbital timescale, 
and Goodman & Oh (1997), in which the efficiency decreases 
almost quadratically. For details on these calculations, see Van 
Kerkwijk et al. (2000). We should stress that at present it is not 
clear that either prescription is reliable; see Goodman & Oh 
(1997) for a discussion. 

The resulting theoretical circularisation timescales are 
listed in Table 3. It is clear that in Goodman & Oh (1997)'s 
prescription the derived timescales are too long by at least an 
order of magnitude, whereas in Zahn (1989)'s prescription a 
timescale is derived that is consistent with the upper limit cited 
above. However, this is only true if the companion is bloated. 
This reinforces the conclusion derived from studies of the cir- 
cularisation of main-sequence binaries that our current theo- 
ries underestimate the efficiency of tidal energy dissipation (see 
Meibom & Mathieu 2004 for a recent summary of the observa- 
tions). 

4.3. The MC-scenario 

The main alternative formation scenario assumes that the sys- 
tem formed through the accretion-induced collapse (AIC) of 
a white dwarf. In this scenario, the companion was filling its 
Roche-lobe at the time of collapse. The kick, that is required to 
bring the system to its present location far outside the clus- 
ter core, increases the orbital separation by only a factor of 
~ 1 .25. Therefore, the companion should still be close to filling 
its Roche-lobe. The only possible way to explain the observed 
faint companion, is that the companion was bloated before col- 
lapse, by irradiation from the accreting white dwarf (as sug- 
gested by Wijers & Paczynski 1993). As it is not clear whether 
this is possible by irradiating only one side of the star (King 
et al. 1996), we regard the AlC-scenario as the less likely for- 
mation model for this system. 

In this scenario, an initial eccentricity of ~ 0.2 is ex- 
pected (Van Kerkwijk et al. 2000). Therefore, the circulariza- 
tion timescale required to explain the current circular orbit is 
roughly the same as for scenario in which the system is formed 
by a triple interaction (^ 2.5 Myr). Even taking into account 
that the companion will be close to Roche-lobe filling, one sees 
from Table 3 that the circularization timescale can be matched 
theoretically only using the prescription of Zahn (1989). Thus, 
our conclusions in §4.2 remain unchanged. 

It is unclear in what precise state the companion is expected 
to be immediately after collapse. For a crude estimate, we as- 
sume that at the moment of collapse, the companion had a ra- 
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dius equal to the Roche limit, and use Eq. 4 to compute the 
critical eccentricity, substituting R c = R L = (l/1.25)/?j n ;. This 
yields, 



e L .. = 0.1383, 



(8) 



where we used the same approximations for the semi-major 
axis and the Roche limit as in the derivation of Eq. 5. This is 
consistent with our measurement of a companion that is bloated 
to a radius smaller than the Roche limit, and thus does not add 
further constraints on which formation scenario is the correct 
one. 

5. Conclusion 

We have used VLT observations to measure the temperature 
and radius of the likely companion of PSR B1718-19. These, 
in conjunction with the companion mass inferred from radius 
and temperature employing pre-main-sequence cooling tracks, 
point to a companion that is bloated, but does not fill its Roche- 
lobe. 

We have shown that a bloated companion is expected in the 
context of the scenario in which the system formed in a a three- 
body interaction in the cluster core less than 10 Myrs ago. With 
the eccentricity that is expected to arise from such an event, 
tidal circularisation can supply the energy needed for a main 
sequence star of the observed mass to expand to the observed 
radius. 

We tried to compare the observed upper limit on the cir- 
cularisation time scale with theory, but were hampered by the 
uncertainty in the extent to which the circularisation efficiency 
is suppressed by the orbital period being shorter than the con- 
vective turnover time. The system constitutes an observational 
constraint on any theoretical account of tidal circularisation. 

The constraints on the properties of the companion could 
be improved by even deeper observations, as these would al- 
low a comparison of the companion with main sequence stars 
of the same color, thus removing the uncertainty in distance 
and reddening. At this point, however, it seems more urgent to 
establish a satisfying theoretical account of the tidal dissipation 
mechanism that is dominant in a fully convective star. 
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